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Abstract 

We study the possibility of relaxing the cosmological bound on the Higgs 
mass coming from the requirement of non-erasure of the baryon asymmetry by 
sphalerons. After reviewing the Standard Model case we obtain this bound in 
two extensions of it: 1) The Standard Model with an additional gauge singlet, 
and 2) The Minimal Supersymmetric Standard Model. Taking fully into account 
all experimental constraints and thermal screening effects we found that the 
situation can be slightly improved with respect to the Standard Model but only 
in case 1) a non negligible region in parameter space exists where the baryon 
washout is avoided and the experimental bounds evaded. 
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Introduction 

The attractive possibility that the observed baryon asymmetry of the universe was 
created at the electroweak phase transition has deserved a lot of of attention in the 
last years [|TJ. It was realized that all the conditions formulated by Sakharov [[| 
(for generating the baryon asymmetry) can be met in the Standard Model (SM). In 
particular, the sphaleron transitions were recognized as unsupressed at high tempera- 
tures and so playing a central role in the mechanisms generating the amount of AB 
we finally observe today 

These sphaleron processes, if unsupressed after the electroweak phase transition, 
will drive the previously created baryon asymmetry to zero, and so its rate V at that 
time 0, U (r ~ exp(—E sp h/T)) should be less than the Hubble expansion rate H to 
go out of equilibrium. Due to the exponential sensitivity of V to the ratio E sp h(T)/T, 
imposing T < H gives the model independent condition ||: 

E sph (T c ) > ^ 



T, 



c 



After relating the sphaleron energy |5| E sp h to the vacuum expectation value of the 
Higgs field at the critical temperature T c (E sp h(T c ) ~ (<p(T c ))/g) the central relation 
Eq. ([I]) translates into 

2 1, (2) 

which tells that the transition should be strongly enough first order so as to avoid the 
erasure of the baryon asymmetry. 

When the dependence of (0(T C )) on the parameters of the Higgs potential is taken 
into account, condition (2) gives an upper bound on the mass of the Higgs boson (in 
the Standard Model where this calculation was originally performed and also in more 
general cases). The study of this bound is the central aim of the work here resumed. 

All along this paper we will say that a phase transition is strongly first order if 
condition @) is satisfied and weakly first order (or even second order) if it is not. 



The Standard Model 

The effective potential at non zero temperature is the central tool in studying the 
electroweak phase transition. In one-loop approximation the temperature-dependent 
potential describes a first order phase transition between the symmetric phase (at high 
T) and the SU(2) x U(l) breaking one. Using a high temperature expansion the form 
of the potential is: 

V{<f>, T) = D(T 2 - T o 2 )0 2 - ET(p 3 + ^A(T)0 4 , (3) 
where D,T Q , E and A(T) are easily calculable. 
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Approximating the critical temperature T c as the temperature Tp at which two 
degenerate minima coexist (T c will be somewhere in the narrow range between Tp and 
the temperature T at which the minimum at = is destabilized) one gets 

(<KTc)) 2E 

T c \(T c y { ] 

We see from this formula that the strenght of the transition is governed by E (which 
gives the cubic term in the potential). Now, E comes from bosonic loops only (from 
n = Matsubara frecuencies) . Each bosonic degree of freedom, with field dependent 
mass Mi((f)) contributes to the potential a term 



T 



3/2 



M/(0) • (5) 



12tt 

In the SM, the main contribution is the one coming from gauge boson loops: 

T 



AV 



12tt 



1 \3/2 n x 3/2' 



4 J V4 



(6) 



Using the value of E obtained from Eq. (^) and relating A to the Higgs mass one gets 
the bound || 

M h < (45 - 50) GeV, (7) 

region which is already ruled out by LEP. 

This bound is obtained within the one-loop approximation to the effective potential 
and as is well known this will only be reliable for <fi > T (the expansion parameter is 
g 2 T 2 1 M 2 (</>)) . For small <fi the perturbative expansion is plagued with infrared diver- 
gences. To cure this, one has to reorganize the loop expansion resumming the most 
important infrared divergent diagrams || (the so called daisy diagrams). The effect 
of this daisy resummation is to shift the mass of the n = modes contributing to the 
cubic term by a thermal mass (Debye screening effect). In the SM this thermal mass 
is of order g 2 T 2 to leading order for the longitudinal gauge bosons, and zero (to this 
order) for the transverse ones. The effect of this thermal mass is to shield the cubic 
term of the longitudinal modes while leaving the transverse modes unaffected and so 
to effectively reduce E by a factor 2/3. The improved bound on the Higgs mass is then 
even lower [PI: 



M h < V2/3(45 - 50) GeV ~ (35 - 40) GeV. (8) 

There exist further refinements of this calculation in the literature. If the screen- 
ing of the transverse modes at subleading order (g 4 T 2 ) is taken into account [10] the 



strenght of the transition is further reduced. On the other hand, some non-perturbative 
effects may give a lower value for the critical temperature ]TTJ translating into a larger 
ratio (0(T c ))/T c . Also there are some two-loop (resummed) calculations |12[ that give 
larger (4>(T C )). 
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Even with these uncertainties the SM has another problems to be able to generate 
the observed baryon asymmetry. In fact the CP violation effects are by far too small 



to give the correct amount of AB and some recent atempts |T3 to solve this problem 



have been shown |L4j] to be flawn. So it seems necessary to move to some extensions of 
the SM to explain the generation of the baryon asymmetry at the electroweak phase 
transition. 

From what we have learned by studying the Standard Model case we can work out 
an estrategy |15[ for the extensions of it to give a strong first order phase transition. 



A look at Eq. (|J) shows that we need a large cubic term in our potential. So, adding 
more bosonic degrees of freedom may help for only bosons contribute to this cubic 
term. The screened mass of such bosons will have the general form 

M 2 (0) = M 2 + A0 2 + n(T), (9) 

where M is an invariant mass, A is a generic coupling of the new bosons to the Higgs 
(here <fi stands generically for the fields driving the transition) and II is the thermal 
mass. For this to give a non-negligible contribution to the cubic term, small values of 
M and II and large values of A are required. Usually, there will exist an experimental 
lower bound for M so that one cannot take it to be arbitrarily small. Concerning the 
coupling A the requirement of perturbativity up to some large scale will set an upper 
bound on it. In a given model one has to take all this conditions into account when 
looking for a region in the parameter space where the phase transition is strong enough. 
Another way of increasing the strenght of the transition has been considered in the 



literature [[15], |16|, [T?| , namely the possibility of decreasing the value of the coupling in 
the denominator of Eq. @) due to one loop corrections (even at zero temperature). 
But usually for this mechanism to give a sizeable effect large couplings and/or many 
bosonic degrees of freedom are needed. 

Before moving to the discussion of extended models some words about the validity 
of the perturbative calculations are needed. Roughly speaking, one can trust the 
(resummed) perturbative calculations in the SM if m\ £ (or equivalently X/g 2 < 1). 
When trying to improve the "sphaleron bound" < mh,criti one can face cases where 
inherit ~ mw rising some doubts about the reliability of the calculation of rrih^rit- in 
the cases we are going to consider the relevant parameter will no longer be X/g 2 . In 
fact g 2 will be substituted by ( 2 in the SM + singlet case and by h 2 in the MSSM . 
As we will be interested in large values of these couplings, for a given Higgs mass the 
expansion parameter will be smaller than in the SM, and the perturbative calculation 
more reliable. 



The Standard Model with a gauge Singlet 

This is the simplest extension of the Standard Model that can in principle improve the 
situation concerning the electroweak baryogenesis problem, as already noted in Ref. 
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The lagrangian of the model is defined as: 



£ = £ 



SM + d^S^S 



M 2 S*S - \ S (S*S) 2 - 2( 2 S*SH*H, 



(10) 



where H is the SM doublet with (H) = <p/ V% (j> ls the classical field, and M 2 , As, ( 2 > 
0, to guarantee that (S) = at all temperatures. 

The temperature dependent effective potential in the one loop approximation was 
studied in Ref. ||15|| . The daisy improvement was performed in Refs. [[18|, [Tj|. The 
potential is a function of the classical field <fi as in the SM, but with a new contribution 
coming from the S bosonic loops of the form (imposing renormalization conditions 
preserving the tree level value of the vev v) 



AV S = g s 



m|(0)T 2 ikf s (0)T m|(0) 



24 



12tt 



64tt 2 



log 



c B T 2 



m|(0) 



where ^5 = 2 is the number of degrees of freedom of the (complex) singlet. mg(0) is 
the field dependent mass of S and M s ((f>) its Debye screened mass. They are given by 



m|(0) 



2 J.2 



Aft 



M 2 + C 



m!(0)+n s (T). 



:i2i 



Here is, to leading order 



n 5 (T) 



(13) 



The thermal polarizations for the rest of SM particles can receive extra contributions 
coming from S loops (see Refs. [|1^, [1(| ) . 

To understand the numerical results presented at the end of this section it is conve- 
nient to study analytically the effective potential assuming as will be the case that the 
bosonic contribution is dominated by the singlet (and this is the reason for the relevant 
expansion parameter to be \/( 2 rather than X/g 2 )- The high temperature expansion 
of the effective potential takes then the form 



V{<f>) = A{T)4> A + B(T)tf + CiT) U z + K 2 (T) 



3/2 



where 



C{T) 



C 3 T 



K 2 (T) 



6/T 

(C 2 + A S )T 2 



3M 2 



3C 2 



(14) 



(15) 



(16) 



From these expressions and the discussion at the end of the previous section is clear 
that the best case for the phase transition to be stongly first order is to have low values 
of M and Xs and large values of ( 2 . This is what was found numerically in Ref. 118]. 
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For large values of M the transition is weaker and eventually, when M 3> T, the 
singlet decouples and the SM result is recovered. The dependence on ( 2 is the most 
interesting one. Imposing the condition that all the couplings in the theory remain 
perturbative up to some high scale A an upper bound for ( 2 at the electroweak scale is 
derived (as a function of course of the scale A and the unknown value of the top yukawa 
coupling ht). This renormalization group study was also performed in Ref. ||18|| . 

Translating the requirement of a strong phase transition (i? S p^(T c )/T c > 45 or 
(0(T C )) > 1.3 with our definition of T c ) into a bound on the Higgs mass M^ lt , one finds 
that in this model one can evade the stringent bound of the SM and go to larger values 
of Mft %t . The numerical value of this critical mass depends on how large the coupling 
( 2 is, or equivalently, on the mass of the top and the scale of new physics A. Assuming 
that this scale is not lower than 10 6 GeV the value of Mfj lt can be of order 80 GeV 
(see Fig. 1) which is higher than the experimental bound. Larger values of this critical 
mass can only be obtained at the price of introducing new physics at lower scales. 



The Minimal Supersymmetric Standard Model 

The Minimal Supersymmetric Standard Model (MSSM) is the physically most mo- 
tivated and phenomenologically most acceptable among the extensions of the SM. 
Concerning the generation of the baryon asymmetry the MSSM allows for extra CP- 
violating phases besides the Kobayashi-Maskawa one, which could help in generating 
the observed baryon asymmetry |19). On the other hand the MSSM has a lot of new de- 
grees of freedom and the nature of the phase transition could be significantly modified 
with respect to the SM. 

The main tool for our study is the one- loop, daisy-improved (for previous studies in 
the one-loop approximation see Refs. 20]) finite-temperature effective potential of 
the MSSM, Kff(0j T). Now the potential is actually a function of two fields: 4>i = R- e ^i 
and 4>2 = Reif^ where Hf and H® are the neutral components of the Higgs doublets 
Hi and H 2 , thus will stand for (0!,0 2 ). Working in the 't Hooft-Landau gauge and 
in the DR-scheme, we can write 



where 



VM,T) = V {<f>) + ^(0,0) + AVi(0,T) + AV daisy (0,T) 



W,0) 
AVi(0,T) 
AV daisy (0,T) 



m\(t)l + ml4>l + 2m30i02 + - +9 (4>l 




m\ (0) 3 
log - - 



Q 2 



y2 



12tt 



^77, [mf(0,T) -777^(0) 



(17) 

(18) 
(19) 

(20) 

(21) 
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The first term, Eq. (|Hj|), is the tree-level potential. The second term, Eq. ([H]), is 
the one- loop contribution at T = 0: Q is the renormalization scale, where we choose 
for definiteness Q 2 = m|, m 2 ((p) is the field-dependent mass of the i th particle, and 
Hi is the corresponding number of degrees of freedom, taken negative for fermions. 
Since Vi(4>, 0) is dominated by top (i) and stop (£1,^2) contributions, only these will 
be included in the following. The third term, Eq. fl2~0|), is the additional one- loop 
contribution due to temperature effects. Here Jj = J + (J_) if the i th particle is a boson 
(fermion), and 

J ± (y 2 ) = J™ dxx 2 log =F e-v 7 ^^ . (22) 

Since the relevant contributions to AVi(0, T) are due to top (£), stops (ti,^) an d 
gauge bosons (W,Z), only these will be considered in the following. Finally, the last 
term, Eq. (pll) , is the correction coming from daisy diagrams. The sum runs over 
bosons only. As usual, the masses m 2 (<p,T) are obtained from the m 2 ((f)) by adding 
the leading T-dependent self-energy contributions, which are proportional to T 2 . 
The relevant degrees of freedom for our calculation are: 

n t = -12, nil =n i2 =6, n w = 6 , n z = 3 , u Wl = 2 , n Zlj =n lL = l. (23) 
The field- dependent top mass is 

m 2 t {4>) = hU 2 2 . (24) 
The entries of the field- dependent stop mass matrix are 

m lM) = «4 + ™?(0) + ^ L (0), (25) 
rr$ s {4>) = rn 2 U3 +m 2 M+Dt R (<j>), (26) 
m 2 x {<f>) = h t {A t <t> 2 + i4 x ) , (27) 

where mg 3 , mu 3 and A t are soft supersymmetry-breaking mass parameters, /x is a 
superpotential Higgs mass term, and 

Di(<f>) = Q-^sin 2 ^)^±^(^-^), (28) 

DIM = f|sin 2 ^^-t^(^-^) (29) 
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are the Z)-term contributions. The field-dependent stop masses are then 



m? (0)+m? (0) 
«\ 2 (0) = — o ± 



\ 



mfj0) -m| fl (0) 



[™x(0)f- (30) 



The corresponding effective T-dependent masses, m~ (0, T), are given by expressions 
identical to (150), apart from the replacement 

-> 0' T ) s + n ^«( T ) • ( 31 ) 



6 



The Ilf R (T) are the leading parts of the T-dependent self-energies of ti,R , 

= A -g 2 J 2 + \g 2 T 2 + ^g' 2 T 2 + l -h 2 t T 2 , 
%(T) = ^ S 2 T 2 + 1,' 2 T 2 + ^ 2 T 2 , 



(32) 
(33) 



where g s is the strong gauge coupling constant. Only loops of gauge bosons, Higgs 
bosons and third generation squarks have been included, implicitly assuming that all 
remaining supersymmetric particles are heavy and decouple. If some of these are also 
light, the plasma masses for the stops will be even larger, further suppressing the 
effects of the associated cubic terms, and therefore weakening the first-order nature of 
the phase transition. 

We expect the stops to play an important role in making the transition stronger as 
they have a large number of degrees of freedom and couple to the Higgses with strength 
h t . On the other hand as they are coloured scalars their thermal mass is of order g 2 s T 2 
as shown above and also they are not protected by any symmetry to have an invariant 
mass. The final importance they have in the strenght of the transition will result from 
the interplay between these opposite properties. 

Finally, the field-dependent gauge boson masses are 



a 2 + a' 2 



2 ' ^ vr/ 2 

and the effective T-dependent masses of the longitudinal gauge bosons are 



(34) 



± 



™2 



m 2 



T) =m 2 w ^) + U WL (T), 
T)= l -\m 2 z {4>) + n WL {T) + Tl BL {T) 



\4 



g 2 - 9' 2 



? 1 + ( p 2 )+u WL (T)-u BL (T) 



+ 



2 



(35) 



(36) 



In eqs. (^3j) and , II vi^ (T) and U Bl (T) are the leading parts of the T-dependent 
self-energies of Wi and B L , given by 



Tl WL {T) = -g 2 T 2 



Yl BL {T) = ^g' 2 T 2 



(37) 



where only loops of Higgs bosons, gauge bosons, Standard Model fermions and third- 
generation squarks have been included. 

We need to analyse the effective potential (|T7|) as a function of <fi and T. Before 
doing this, however, we trade the parameters m 2 ,m|,m 2 appearing in the tree-level 
potential fll8|) for more convenient parameters. To this purpose, we first minimize the 
zero-temperature effective potential, i.e. we impose the vanishing of the first deriva- 
tives of Vo(4>) + Vi(0, 0) at (0i,02) = {vijVz), where (t>i,t>2) are the one-loop vacuum 
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expectation values at T = 0. This allows us to eliminate m\ and m\ in favour of 
m| and tan/3 = v 2 /vi in the standard way. Moreover, m| can be traded for the one- 
loop-corrected mass m\ of the CP-odd neutral Higgs boson [51]. Therefore the whole 



effective potential ( |17|) is completely determined, in our approximation, by the param- 
eters (m^, tan/3) of the Higgs sector, and by the parameters (m t , rriQ 3 , mjj 3 , p, A t ) of 
the top/stop sector. The same set of parameters also determines the one- loop-corrected 
masses and couplings of the MSSM Higgs bosons. 

The next steps are the computation of the critical temperature and of the location 
of the minimum of the effective potential at the critical temperature. We define here 
T as the temperature at which the determinant of the second derivatives of V e s(4>, T) 
at <t> = vanishes: 



det 



d 2 V^U,T, 



= . (3£ 

01,2=0 



It is straightforward to compute the derivatives in eq. (^) from the previous formulae; 
the explicit expressions are given in Ref. p^| . 

Once eq. (|38|) is solved (numerically) and To is found, one can minimize (numer- 
ically) the potential V e s(4>, T ) and find the minimum [t>i(T ), t>2(T )]. The quantity 
of interest is indeed, as will be discussed later, the ratio v(T )/T , where v(T ) = 

y/vKT )+v&T Q ). 

Before presenting the numerical results we discuss the experimental constraints on 
the parameters of the top/stop sector and of the Higgs sector. We treat mg 3 , m U3 
and the other soft mass terms as independent parameters, even if they can be related 
in specific SUGRA models. We want to be as general as possible and in this way the 
region of parameters we are able to exclude will be excluded in any particular model. 

Direct and indirect searches at LEP imply |23] that m~ b > 45 GeV, which in turn 



translates into a bound in the (m<g 3 , tan/3) plane. Electroweak precision measurements 
[24] put stringent constraints on a light stop-sbottom sector: in first approximation, 
and taking into account possible effects f25j of other light particles of the MSSM, we 



conservatively summarize the constraints by Ap(t, b) + Ap(t, b) < 0.01. 

We finally need to consider the constraints coming from LEP searches for super- 
symmetric Higgs bosons |5J|. Experimentalists put limits on the processes Z — > hZ* 



and Z — > hA, where h is the lighter neutral CP-even boson. We need to translate 
these limits into exclusion contours in the (mA,tan/3) plane, for given values of the 
top/stop parameters. In order to do this, we identify the value of BRiZ — > hZ*), 
which corresponds to the limit > 63.5 GeV on the SM Higgs, and the value of 
BRiZ — > hA), which best fits the published limits for the representative parameter 
choice m t = 140 GeV, rriQ a = mu 3 = m = 1 TeV, A t = p = 0. We then compare those 
values of BRiZ — > hZ*) and BRiZ — > hA) with the theoretical predictions of the 
MSSM, for any desired parameter choice and after including the radiative corrections 



associated to top/stop loops ^6], ^TJ. Of course, this procedure is not entirely correct, 
since it ignores the variations of the efficiencies with the Higgs masses and branching 
ratios, as well as the possible presence of candidate events at some mass values, but it 
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is adequate for our purposes. 

We now present our numerical results, based on the effective potential of eq. (JTj 
concerning the strength of the electroweak phase transition and the condition for pre- 
serving the baryon asymmetry. 

Particularizing to the MSSM the studies of sphalerons in general two-Higgs models 



[ f2~7jj , we obtain that 

-'sph V- 1 ) — sph 

where, in our conventions, 



Kh SM (T) < E™{T) , (39) 



ggg) 4V27T „j\ eS (T)\v(T) 

and B is a smoothly varying function whose values can be found in Ref. ||. For 
example, ^lO" 2 ) = 1.67, ^(Kr 1 ) = 1.83, 5(1) = 2.10. It can also be shown that 

^1<^1<^>, (41) 

+ D J-C J 

where T c is the actual temperature at which the phase transition occurs, satisfying the 
inequalities 

T < T c < T D , (42) 

if T is defined by ( j38|) and To is the temperature at which there are two degenerate 
minima. 

Finally, the corrections in E~fi due to g' ^ have been estimated and shown to be 
small ||28|| . Therefore, a conservative bound to be imposed is 

a. y^ffg^ i. (43, 

T 45/7 

The last point to be discussed is the determination of the value of A e fr(T ) to be 
plugged into eq. (fUjf). The 5-function we use is taken from Ref. where the sphaleron 
energy was computed using the zero-temperature 'Mexican-hat' potential, V = j((f> 2 — 
v 2 ) 2 . The sphaleron energy at finite temperature was computed in Ref. ||29|| , where it 
was proven that it scales like v(T), i.e. 

E^ h (T) = E^(0) V -^l, (44) 

with great accuracy. Therefore, to determine the value of A e fr(T ) we have fitted 
V e g((j), Tq), in the direction of the minimum, to the one-dimensional approximate po- 
tential, 

\4ff(0,T o ) ~ ^A cfr (T o )[0 2 - v 2 (T )} 2 , (45) 
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The value of X e g obtained from ((45]), 



Aeff(To) = 4 — — , (46) 

where all quantities on the right-hand side are calculated numerically from the potential 
of eq. flT?D, is then plugged into eq. ( f43|) to obtain our bounds. Of course, the quality 
of the fit is good only for values of <p ~ v (To) but this is precisely the region of interest 
to determine the sphaleron energy. 

Our numerical results are summarized in fig. 2, in the (tua, tan/3) plane and for two 
representative values of the top quark mass: m t = 130 GeV (fig. 2a) and m t = 170 GeV 
(fig. 2b). In each case, the values of the remaining free parameters have been chosen 
in order to maximize the strength of the phase transition, given the experimental 
constraints on the top-stop sector. Notice that arbitrarily small values of mu 3 cannot be 
excluded on general grounds, even if they are disfavoured by model calculations. Also, 
we have explicitly checked that, as in Ref. [^D[, mixing effects in the stop mass matrix 



always worsen the case. In fig. 2, solid lines correspond to contours of constant R: one 
can see that the requirement of large values of R favours small tan/5 and rriA ^> ™>z- 
The thick solid line corresponds to the limits coming from Higgs searches at LEP: for 
our parameter choices, the allowed regions correspond to large tan /3 and/or tua 3> mz- 
For reference, contours of constant (in GeV) have also been plotted as dashed 
lines. One can see that, even for third-generation squarks as light as allowed by all 
phenomenological constraints, only a very small globally allowed region can exist in the 
(iriA, tan (3) plane, and that the most favourable situation is the one already discussed 
in Ref. |31J. More precisely, the region that is still marginally allowed corresponds to 
ttia 3> m z , tan/3 ~ 2, stop and sbottom sectors as light as otherwise allowed, a heavy 
top, and a light Higgs boson with SM-like properties and mass nth ~ 65 GeV, just 
above the present experimental limit. A less conservative interpretation of the limits 
from precision measurements, the inclusion of some theoretically motivated constraints 
on the model parameters, or a few GeV improvement in the SM Higgs mass limit, would 
each be enough to fully exclude electroweak baryogenesis in the MSSM. 
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Figure captions 

Fig.l: Plot of (0(T c ))/T c versus m h for \ s = 0, M = 0, m t = 90, 150 and 175 GeV and 
( corresponding to A = 10 6 GeV. 

Fig. 2: Contours of R in the (m A ,tan(3) plane, for the parameter choices: a) m t = 130 
GeV, m Q3 = 50 GeV, m Us = {m~ t ~ 130 GeV, m~ H ~ 50 GeV), A t = fx = 

0; b) m t = 170 GeV, m Qs = 280 GeV, m Us = (m t ~/~ 330 GeV, m~ tR ~ 170 
GeV, mi ~ 280 GeV), A t — fi — 0. The region excluded by Higgs searches at 

LEP is delimited by the thick solid line. For reference, contours of constant 

(in GeV) are also represented as dashed lines. 
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